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Hamilton- Jacobi theory for general relativity provides an elegant covariant formulation of the 
gravitational field. A general 'coordinate-free' method of integrating the functional Hamilton- Jacobi 
equation for gravity and matter is described. This series approximation method represents a large 
generalization of the spatial gradient expansion that had been employed earlier. Additional solutions 
may be constructed using a nonlinear superposition principle. This formalism may be applied to 
problems in cosmology. 

I. INTRODUCTION 

General relativity was formulated to describe the gravitational field in a manner which was independent of any 
particular choice of reference frame. It is invariant under 'general coordinate transformations'. However, when one 
actually solves the Einstein field equations, say in a cosmological setting, one typically invokes a particular choice of 
coordinates. As an improvement to this situation, one would prefer to utilize a 'coordinate-free' method of computing 
solutions to Einstein gravity. In fact, Hamilton-Jacobi (HJ) theory for general relativity provides such a description 
of the gravitational field Q , || . In solving the HJ equation, one need not specify the choice of time hypersurface 
nor the spatial coordinates. 

Coordinate- free descriptions have proven useful in many fields of theoretical physics. Vector analysis for Euclidean 
space is a good example. Although one normally associates coordinates with vectors, it is possible to interpret a 
vector geometrically in a coordinate free way using a magnitude and a direction. In addition, by using his bra and 
ket notation, Dirac was able to formulate quantum mechanics in a form which was independent of the choice of 
basis states. The bra could be interpreted as an abstract vector without referring to a particular basis. 

It has been known since the early 1960's that the HJ equation for general relativity does not refer to the lapse 
and shift parameters which characterize the gauge freedom of gravity j^]. Hence the HJ equation is the natural 
starting point for a coordinate-free analysis of gravity. However, the HJ equation is a nonlinear, functional differential 
equation which characterizes an ensemble of universes, superspace. It was generally believed that superspace was 
too complicated to solve in its entirety. Beginning with Misner Q, researchers were content to solve homogeneous 
minisuperspace models where one considered only a finite number of degrees of freedom for the gravitational field. 
These investigations have been proven to be quite fruitful in quantum cosmology (see, e.g., Louko and Ruback [Q). 
(For a discussion of future trends in quantum cosmology, consult Hartle Q as well as Barvinsky and Kamenshchik 
H.) However, there are some questions that these models cannot address since they do not include inhomogeneities. 
For such models, the choice of time hypersurface is degenerate: a hypersurface of uniform <j> is the same as that of 
uniform scalar factor a. In order to obtain a deeper understanding for the nature of time in a semiclassical setting flTo| , 
|pr| , it is essential to include the role of inhomogeneities. (However, there are many additional obstacles in addressing 
the question of time in a quantum context ]T^| , Ji"3| |, H].) 

It is somewhat surprising that one can actually obtain solutions to semiclassical superspace using a series approxi- 
mation method. The prototype solution for the Hamilton-Jacobi (HJ)equation for general relativity utilized a spatial 
gradient approximation. It was developed by Salopek and Stewart (TJ and advanced further by Parry, Salopek and 
Stewart |2|. In effect, they demonstrated that one can decompose superspace into a discrete sum of minisuperspaces. 
Here, a general method will be given for constructing (HJ) solutions which utilize various other series approximations. 
These techniques may be applied profitably to cosmology and other areas of astrophysics. 

A Hamilton-Jacobi description of general relativity is attractive for several other reasons: 

(1) Primitive formulation of quantum theory for the gravitational field. In the semiclassical approximation, the 
wavefunctional is given by a phase factor, \& ~ e lS l h , where Planck's constant h~ is assumed to be tiny. The generating 
functional S then satisfies the HJ equation. If S is real, then one is describing classical phenomena. If S is complex, 
one may describe quantum phenomena such as tunnelling or the initial wavefunction of the universe. It is generally 
believed that the fluctuations for structure formation as well as microwave background fluctuations are generated 
during an inflationary epoch. For such models, it is absolutely essential to quantize the gravitational field, including 
both scalar and tensor modes (see, for example, The current status of inflation after the detection of the 
cosmic microwave background anisotropy is discussed in refs. Jl6[ . 

(2) Systematic solution of constraint equations. In superspace, the momentum constraint may be given a simple 
geometric interpretation and it is easy to solve. The energy constraint fully characterizes the dynamics of the gravi- 
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tational system, and it is difficult to integrate directly. As will be demonstrated in the present work, it will be solved 
using a series approximation. 

(3) Nonlinear analysis of gravity. In the quasi-nonlinear regime, one may employ HJ theory to analyze the formation 
of cosmological pancakes ]17j during the matter-dominated epoch of the Universe [Q , jl9| . 

Spatial gradient expansions in general relativity have a long history p3]. Within a HJ context, a spatial gradient 
expansion has been applied successfully to cosmology including a detailed computation of microwave background 
fluctuations and galaxy-galaxy correlations function arising from inflation , |2^] . Soda et al [^3| have generalized the 
expansion to encompass Brans-Dicke gravity, and Chiba |2J] has formulated the gradient expansion for n-dimensional 
Einstein gravity. The HJ equation for long- wavelength fields has often been invoked in an attempt to recover the 
inflation potential from cosmological observations |25[ , [^6| , (2^] . 

In Sec. II, the HJ equation for a scalar field interacting with gravity is presented along with the analogous equation 
for a dust field interacting with gravity. I quickly review the spatial gradient expansion of the generating functional 
in a way which is easy to generalize to other situations. In Sec. Ill, the HJ equation including a scalar field is solved 
using a Taylor series in the scalar field <fi. In Sec. IV, the analogous method is applied for a dust field, which describes 
collisionless, pressureless matter. In Sec. V, the Superposition Principle for Hamilton- J acobi theory allows one to 
construct complicated solutions of the HJ from known solutions which depend on a parameter. Conclusions and a 
summary follow in Sec. VI. 

(Units are chosen so that c = 8irG = 8ir /m 2 P = 1. The sign conventions of Misner, Thorne and Wheeler |2f| will 
be adopted throughout.) 



II. HAMILTON- J ACOBI EQUATION FOR GENERAL RELATIVITY WITH MATTER 



In the present work, two situations will be considered: (1) a scalar field interacting with gravity and (2) a dust field 
interacting with gravity. 



A. Scalar field interacting with gravity 

The action for a single scalar field interacting with Einstein gravity is 

I = J afx^—g (IWR-lgTd^^-V^ ; (2.1) 
( 4 ) R is the Ricci scalar of the space-time metric g^ v . For simplicity, the scalar field potential is assumed to be 

V(<t>) = \m 2 4> 2 + V Q (2.2) 

which describes a massive scalar field with a cosmological constant term. (In general the solution methods described 
in this paper will work for all potentials which are regular at = 0.) 
In the ADM formalism the line element is written as 

ds 2 = (-N 2 + j ij NiNj) alt 2 + 2N l0 lt dx l + j^dx* dx j , (2.3) 

where TV and Ni are the lapse and shift functions respectively, and 7^ is the 3-metric. (For a elegant review, consult 
D'Eath [ESI). One can then rewrite the action in Hamiltonian form, 



d 4 x 



(n<t><p + Tr' lj % — NH — N l n^ . (2.4) 



where TC and Hi are the energy and momentum densities, respectively. In the Hamilton- Jacobi formalism, one replaces 
the momenta by functional derivatives of the generating functional, S["f a b(x), </>(x)], 

<hfij(x) 54>{x) 

the generating functional associates a complex number for each field configuration <p(x) on a space-like hypersurface 
whose geometry is described by the 3-metric jij(x). The energy and momentum constraints are now given by: 
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= H(x)=j-^ 2 [2j ik (x)j j i(x)-j ij (x)j k i(x)] 



and 



Sjij (x) S^/ki {x) 

b~"* (im) 2 + ~< 1I2V ^')\ - ^" 2 R + ^"VVA , (2.6a) 



" = >" ' = ' , + J&p«M + <2 ' 6b> 



The Hamilton-Jacobi equation (2.6a) governs the evolution of the generating functional, S = S[j a b{x),(j)(x)], in 
superspace. It is quadratic in the momenta of fields. The momentum constraint is linear in momenta. Higgs j^] 
showed that it legislates that the generating functional is invariant under reparametrizations of the spatial coordinates 
('spatial gauge-invariance'). It may be solved, for example, by assuming that the generating functional is an integral 
of some function of the curvature, say J cPxj 1 / 2 f(R), or some other combination of spatial invariants. 



B. Dust field interacting with gravity 



The case of dust interacting with gravity is of high interest in cosmology. It is generally believed that most of 
the Universe consists of dark-matter whose dynamics may be described by a dust field which is pressureless and 
collisionless. Secondly, hypersurfaces of uniform \ in general relativity come closest to describing Lorentz frames that 
proved so useful in flat spacetime. In fact, a Lorentz frame may be considered to be a properly synchronized collection 
of dust particles. It may be necessary to introduce dust to interpret a quantum theory of the gravitational field 
@. 

The action 



fJ 
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WR--n(g'»'d M x&X+l) 



(2.7) 



for a dust field, x> interacting with gravity, is similar to that of a scalar field. The new ingredient is the rest number 
density n = n(t,x) which is a Lagrange multiplier that ensures that the 4-velocity 



satisfies U^Un 



u» = -sTx,u 

T. Hence \ may be interpreted as a velocity potential. In Hamiltonian form, the action is 



d 4 x (tt x x + n ijA fij -NU- N'U,) . 

where 7r x is the canonical momentum of the dust field. Replacing the momenta by functional derivatives of S, 

SS „, , SS 



n %J (x) = 



TT X (x) 



§x(x) 



(2.8) 



(2.9) 



(2.10) 



the energy and momentum constraint equations become: 

= H(x)= 7 -^-^ 



= 7U(x) 



5jij(x) Sju(x) 
r " SS 

ss 

tk $7kj (x) 



[2ja(xh jk {x) - 7«i(ar)7fcj( a; )] 



S X (x) 2 
SS 



S^fki(x) 



Ikl,: 



SS 

Sx(x) 



X,: 



(2.11a) 
(2.11b) 



In contrast to the Hamiltonian for a scalar field (2.6a), the energy constraint (2.11a) is linear in the canonical 
momentum of the dust field. In a quantum context, the Hamiltonian constraint for dust and gravity is thus very 
similar to the Schrodinger equation (or its covariant generalization, the Tomonoga-Schwinger equation |3C| ) that has 
been so successful in flat space-time. 
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C. Review of Spatial Gradient Expansion 



The spatial gradient approximation for the HJ equation (2.6a) with a scalar held will be quickly reviewed. The 
essential aspects will highlighted in order to illustrate what must be done in other situations. 
One expands the generating functional 

oo 

S[lab(x), ^{x)] = <s( 2 ") (spatial gradient expansion) (2-12) 

in a series of terms according to the number of spatial gradient terms that they contain. It is quite important that 
each term in the series expansion satisfy the momentum constraint, 

o> > / SS*- 2 ^ \ (SS( 2n ) (5S( 2 ™) 

o = n? n \x) = -2 httjP-pr) +rn^ + Srr^- (2.13) 

\ ojkj{x) J . O'fki(x) d(p(x) 

The zeroth order term 

5 (0) = -2 /" d 3 x^ 2 H((f>) , (2.14) 
is the simplest such term that one can imagine; it contains no spatial gradients where the function H((f>) satisfies 

H 2 = l( d -^\ +\v{<f>). (2.15) 



3 V d<j> J 3 

The volume element cfixj 1 ^ 2 appearing in eq.(2.14) is obviously invariant under spatial coordinate transformations. 
The second order term contains two spatial gradients and is an integral over the 3-curvature and a term quadratic in 
spatial derivatives of <j>: 

S (2) = f 1 l/2 d 3 x ^ J[m + R{(j)) ^b^^ . (2 16) 



J and K are arbitrary functions of 4> which are chosen to satisfy the second order HJ equation. The higher order 
terms proceed along similar lines: e.g., S 1 - 4 ^ consists of all invariant terms with four spatial derivatives. 

What precisely is the expansion parameter in the gradient expansion? To what does the index 2n refer? The index 
2n is related to the conformal weight of the functional S^ 2n \ To clarify this point, it is useful to introduce a scaling 
factor, s. If one rescales the 3-metric using the homogeneous conformal factor, s, 

lab{x) -> S 2 Jab(x) (2-17) 

one finds that 

5 ( 2 ") [s 2 lab , 0] = S < 3 - 2 ") [ 7o6 , 0] . (2.18) 

Hence the gradient expansion is an expansion of the generating functional in powers of the scaling factor s: 

00 

S[s 2 lab , 0] = s (3 - 2n) <S (2n) [lab, 0] , (2.19) 

where 6>( 2 ™) [j a b, are simply the coefficients of s^ 2n \ At the end of the calculation, one sets the scaling factor to 
unity because it is just a counting parameter. 

The interpretation of the spatial gradient expansion as an expansion in powers of the conformal weight may be 
trivially extended to other situations. For example, instead of rescaling the metric, one may choose to rescale the 
matter field, <j>(x) — » s<j>(x) (or x( x ) ~ * s xi x )) ancl then expand in powers of s. Such a simple adjustment leads to 
radically different forms of the generating functional as will be demonstrated in Sec. Ill and Sec. IV. (In order to 
reduce the introduction of extraneous notation, one in practice expands in powers of the desired field, and typically 
foregoes any mention of the scaling factor s.) 
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III. TAYLOR SERIES EXPANSION IN THE SCALAR FIELD 



For the case of a scalar field interacting with gravity, we will now consider a solution for the generating functional 
in the form of a power series of the scalar field: 

oo 

S[-y ab (x),cj>(x)} = S(m) ( series in C 3 ' 1 ) 

m—0 

The zeroth order term will be assumed to be independent of <p, but otherwise, it is an arbitrary functional of the 
3-metric, 

S<°> =S^[ lab {x)] , (3.2a) 
which is invariant under reparametrizations of the spatial coordinates: 

/ 8S^ \ SSW 

= -2 h ikJ pr + T— T-xTfcW ■ 3 ' 2b 



One substitutes the series into the HJ equation 



A. Equations for Scalar Field 

and collects terms of like order, to find: 



2.6a 



6SW 
J<j>(x) 

ssm ss (2 



8(j)(x) S(j)(x) 
6SW 5S& 



jR - 27 Vq + -2 [2ju(x)jjk{x) - Jijix^kiix)] - — pr- p- , (zeroth order terms) , 

7^j \X) o'Ykl \X ) 

-2 [2rfu(x)'yjk(x) - Jij (x)7jm {x)\ : _ f ^ : _ ,^ , (first order terms) 



5"fij(x) 8jki(x) 



<5<Hz) Jc/>(x) 2 V<^0) 



— TO 

2 
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(second order terms) 



[2jii{x)^ jk {x) - lij[x)^ki{x)\ I t —~ — + 2- 



It is straightforward to derive a general expression for higher order terms in analogy with ref. 



(3.3a) 
(3.3b) 
(3.3c) 



B. Solutions for a Scalar Field 

The generating functional for any order may be written in terms of a recursion relation. The first few terms 



S W=sW[ lab (x)} 



5(1) 
5(2) 

5 (3) 



d s x^' 2 4>{x) 



R-2V Q -2^ 1 [2j u (x)j jk (x) - %j(x)^ kl (x)] 



8SW «js(o) 



d 3 x(f>(x) [2j u (x)7jk(x) - li]{x)lki{x)] 



8S<® 6SW 1 8SW 



1/2 



Sjij(x) 8^m(x) / 6<f>(x) 



d 3 x<i)(x) 



1 (8S^ 



2 \8<j){x) 



7 2 1 2 7 aft / / 



[2jii{x)^j k {x) - 7yW7fe«WJ I ~ — ~ — + 2 " 



(57y (x) (57fez (x) Sjij (x) 8j ki (x) 



8<f>{x) 



(3.4a) 
(3.4b) 

(3.4c) 
(3.4d) 



The validity of these these formulae will be demonstrated by deriving iS^ 2 ). The functional equation ( |3. 3b ) defining 
iS^ 2 ) may be rewritten in the form 
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5<p(x) 



nr. , s , « , , ( m ^ (0) MM ,8SW 

o^ij[x) • 0"{ k i{x) 64>{x) 



(3.5) 



Given the right hand side is known. Hence eq.( |3.5| ) has the form of an infinite dimensional gradient which may 
be integrated using a contour integral || in cf> field-space (the 3-metric is held fixed during such an integration). For 
the same reasons given in earlier work [ ^Cj[ , one may choose an arbitrary contour of integration, and I will use a 
straight line parameterized by 



(j>(x) = t<f>(x) , < t < 1 . 
to connect the origin, 4>o{x) = to the 'final point' (f>f(x) = 4>{x): 



= / d 3 x t dtd>{x)^- 
Jo 5cj)(x) 



(3.6) 



(3.7) 



Aft er c ounting powers of t, the integral over t gives 1/2, and hence is given by e q.(|3.4cp . Integrability of in 
eq^^) is guaranteed if the previous order terms, iS^ ) and are gauge-invariant fiof ^ 



1. First Example for a Scalar Field 

We will now give some explicit examples. Since S^ 3 ' [y a b(x)] is arbitrary, the possibilities are limitless, but a sampling 
is instructive. 



S<°> = -2H J d 3 x7 1 / 2 , H = VW, 
SW = [ d^x^^R 1 ' 2 , 



S<V = -H J d'xj 1 / 2 



(3.8a) 
(3.8b) 

(3.8c) 



No spatial derivatives of <fi appear in . However, they do appear in , and they cannot be removed by integration 
by parts. 

Unlike the spatial gradient expansion, the above series is not analytic since contains a square root of R. If 
R < 0, S is imaginary, which is classically forbidden. If the sign of the square root is chosen accordingly, these 
solutions would be exponentially suppressed in a quantum analysis: 



exp 



-2 J d 3 xj^ 2 <py/=R 



, (for </>>()). 



(3.9) 



2. Second Example for a Scalar Field 



Another simple example arises if we take S(°> = 0: 
=0, 

s« = / dW /2 0(^-2y o ) 1/2 , 



S<»=Q, 



where 



y/R ~ W 

5SQ-) 7V2 



6jij(x) 5~fki{x)_ 



(R - 2Voh ab - R ab + D]ab ~ J ab D 2 } U (R ~ 2Vo) 



-1/2 



5lab{x) 2 

If R > 2Vo, then describes the classically forbidden regime. 



(3.10a) 
(3.10b) 
(3.10c) 
, (3.10d) 

(3.10e) 
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IV. TAYLOR SERIES EXPANSION IN DUST FIELD 



A. Equations for Dust Field 

I will now consider an expansion of S in powers of 

oo 

S[jab{x),X( x )} = ^2 S(m) ( Seri(3S U1 X)- 



(4.1) 



This expansion is quite general, and it will be applicable in most instances except when one encounters a singular 
point (which will require a separate treatment). The zeroth order term will be assumed to be independent of x : 



S< >=S<°> [7.6(a)]. 
Expanding all terms in powers of x, one obtains the following equations: 

1 I/a „ _,/,,„ , * / \ , „ ss^ ssw 

Mx) = 2 7 R ~ 7 [2la {x)ljk {x) ~ 7y (x)7fei ' 

= -27 7 L27iiW7jfe(ar)-7«W7iwWj 



#x(x) 
<JS< 3 > 



57y(a) ^(x) 



SS W 1 / la\ ^ (iJ -1/2 r„ / x / n 

= "2 (XloX 1 J - 7 F^uWTjkW - lij{xhki{x)} 



55W <5S« <fs(°) 
2- 



5 7ii (x) <57fej (x) (x) <5 7fc z (x) 



(4.2) 

(4.3a) 
(4.3b) 
(4.3c) 



B. Solutions for Dust Field 



The above equations (4.3) may be integrated immediately 
5(°) = 5(°)[ 7 a6(x)], 
= [ d 3 x X (x) 



2 ojij(x) 07fe;(x) 



5 (2) 



d 3 x x(x) 7 1/2 [27,7 (x)-fjk (x) - 7ij (x)7 fe( (x)] 



6SW SSW 



8SW 

Sx(x) 

- 7 "V2 [2 7a (a;)7 ifc (x) - 7 4i (x)7 fez (x)] 



<5«S« SS^ 55(a) 
2- 



57ij (x) c>7fci (x) 57y (x) <$7 M (x) 



(4.4a) 
(4.4b) 

(4.4c) 
(4.4d) 



One may interpret the higher order terms, S^ 1 ' , S^ 2 ' , . . ., as describing the evolution in time x(x) of the the initial 
state S(°)[7 a6 (x)]. 



= C / d 3 x 7 1/2 ■ 



1. First Example for a Dust Field 
If <S(°) is proportional to the volume of a given 3-geometry, the first few terms are: 

3 x 7 1 

S« = /dW /2 X [f + 

2 /i? 9C 2 \ 1 , fc 



S< 2 ) - C / d 3 x 7 1 / 2 



(4.5a) 
(4.5b) 
(4.5c) 
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2. Second Example for a Dust Field 



As a second example, consider a series whose first term is proportional to an integral of the spatial curvature: 

5 (0) =E I d 3 x^ 2 R, (4.6a) 



S«= / dV /2 x 



i? 



- 2E Z [ R ab R ab ~^R' 



(4.6b) 



V. SUPERPOSITION OF HAMILTON- J ACOBI SOLUTIONS 



One of the very attractive features of quantum mechanics is the principle of linear superposition: if ip\ and 02 are 
solutions of the quantum theory, then ip(y) = ipi + 02 is a solution as well. This principle has no direct classical 
interpretation. Nonetheless, there are situations where one can indeed construct additional solutions to the HJ 
equation from other known solutions. One can enunciate a Principle of Superposition for Hamilton- J acobi theory 
which is motivated by a semiclassical treatment of the quantum theory. 

For a quantum system with configuration variable, y, suppose that one is fortunate enough to find a solution to 
the Schrodinger equation, vb(y\a) which depends on a continuous parameter a. Then any linear superposition of these 
solutions is also a solution: 

0(2/) = J ' daw(a)ifj(y\a) , (5.1) 

where the weighting function w(a) is arbitrary. If we work in the semiclassical limit, fi — > 0, then rp(y\a) and w(a) 
may be approximated by phase factors, 



■0(2/1 a) 



JS(y\a)/h 



w(a) 



— P ig{a)/h . 



(5.2) 



S(y\a) is then a solution of the HJ equation which depends on a parameter a. The resulting integral ( p.l[ ) may 
approximated using the stationary phase approximation, 

V>(y) - exp[i(S(y\a) + g{a))/a] (5.3a) 
where a = a(y) is now chosen so that the phase of the integrand has a maximum or minimum, 

= ^-[S(y\a)+g(a)}. (5.3b) 
oa 

Hence the principle of linear superposition in quantum mechanics and the stationary phase approximation lead to 
the Principle of Superposition for Hamilton- Jacobi theory: 

If S(y\a) is a solution of the HJ equation which depends on a continuous parameter a, then 

T(y) = S(y\a) + g(a) (5.4a) 
is also a solution provided that a = a(y) is determined by the stationary phase condition, 

dS(y\a) dg(a) 

= — 1- — — . (5.4b) 

oa oa 

The proof given above was motivated by the quantum theory. In a gravitational context, it is not clear that a 
consistent quantum theory exists (at present). However, one can verify the principle totally within a classical context 
by noting that 



dT{y) _ dS(y\a) 
dy dy 



dS(y\a) da dg(a) da 
da dy da dy 



(5.5) 



The last two terms vanish by virtue of the stationary phase condition (5.4b), and thus a derivative of S(y) with 
respect to y coincides with a derivative of S(y\a) with respect to y (holding a fixed): 



dT(y) _ dS(y\a) 
dy dy 



(5.6) 



Since S appears in the HJ equation only in terms of its derivatives with respect to the configuration variables, the 
principle is justified. 
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A. Applying the Superposition to the HJ Equation for Gravity 



The superposition principle leads to some rather exotic solutions of the HJ equation for general relativity. To 
illustrate the basic principles, we will consider a massless scalar field, m — with vanishing cosmological constant, 
Vq = 0, which has the following spatial gradient expansion solution for the HJ eq.(2.6a), 



S[ lab (x),<f>(x)]=sW +5 (2) +..., 
of which the first two terms in a spatial gradient expansion are 

S<°> =-2C J A 7 1/2 e^, 

S (2) = J fxfW \L e ~^ (iRt + \r b M^ +ERf 



(5.7a) 

(5.7b) 
(5.7c) 



C and E are homogeneous but arbitrary parameters. The new 3-metric j a b is related to 7 a b by a conformal transfor- 
mation: 



fab = e 



2<f>/V6, 



"fab ■ 



Rf is the Ricci Scalar associated with / a f,. 

Consider a new solution T of the HJ equation given by the sum 

T = S-\C, 

where T satisfies the stationary phase condition with respect to C: 
Using a spatial gradient expansion, one can solve for C: 



C = C -±J dVVa e -4*/VS (l R f + If*^^ 



+ 



where the functional Co is given by 



Co 



2 / d 3 a; 7 1 /2 e V3/2*. 



T admits the following spatial gradient expansion: 

T = T (0) + T (2) 

with 



T< 2 > =E J d 3 xf/ 2 Rf + J d 3 x f 1 ' 2 e" 4 ^ \-^ Rf ~ I (J cPx^l^e^T^ 



(5.8) 



(5.9a) 



(5.9b) 



(5.10a) 



(5.10b) 



(5.11a) 



(5.11b) 



(5.11c) 



One may verify directly that eq.( 5.11a ) is solution of th e HJ cq.( 2.6a ). 

Unlike the original gradient expansion solution ( |5. 7a ) for S, the new solution T is no longer an integral over local 
terms, but contains products and quotients of local integrals. Hence very complicated solutions of the HJ equations 
may be generated using the nonlinear superposition of local solutions. 
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VI. CONCLUSIONS 



A coordinate- free analysis of general relativity possesses distinct advantages over traditional methods of solving Ein- 
stein's equations. When solving the field equations, one typically makes arbitrary gauge choices. In many situations, 
the optimal choice of gauge is not very clear, and a poor choice of gauge can complicate the analysis significantly. 
In the coordinate-free method expounded here, one may solve the functional HJ equation without making any gauge 
choices. However, applications of the HJ equation to physical problems have previously been hampered by the lack of 
mathematical tools. Some general and useful techniques have been developed which will be applied later in cosmology. 
However, the functional approach is a radically different way of describing gravity. More work is required in order to 
develop a more intuitive understanding of the method. 

Although one can obtain exact general solutions for gravitational superspace in two spacetime dimensions (see, for 
example, |3^] ), it is doubtful that one may construct such solutions in four spacetime dimensions. One must resort to 
some approximation method. The spatial gradient expansion was the prototype solution of the H J equation for gravity 
and matter. It is basically a Taylor series expansion in the conformal weight factor of the 3-metric. By expanding 
the generating functional in terms of other fields, one may construct numerous other solutions to the HJ equation for 
general relativity. Explicit expansions in either a scalar field </> or a dust field \ were demonstrated explicitly. These 
solutions describe the evolution of some gauge-invariant initial state [y a b(x)] as a functional of the matter field. 
Integrability of these solutions is guaranteed by spatial gauge-invariance. 

Many of these solutions depend on continuous parameters. One can construct additional solutions by using the 
Superposition Principle for Hamilton- J acobi Theory which is motivated by the superposition principle in quantum 
mechanics in conjunction with the stationary phase approximation. One can in effect construct complicated solutions 
by integrating over the continuous parameters. 
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